Fully coupled hydraulic fracture simulation using the improved element partition method     

Shujun Peng  Zhennan Zhang  Jianye Mou  Bing Zhao  Zhiyuan Liu  Jiading Wang 《国际地质力学数值与分析法杂志》2019,43(1):441-460

The improved element partition method (IEPM) is a newly developed fracture simulation approach. IEPM allows a fracture to run across an element without introducing extra degrees of freedom. It can also simulate any number of fractures in a prescribed mesh without remeshing. In this study, the IEPM is extended to hydraulic fracture simulation. First, the seepage and volumetric storage matrix of a cracked element are derived using virtual nodes (the intersection points of a crack with element edges). Subsequently, the fully coupled hydromechanical equation is derived for this cracked element. To eliminate the extra degrees of freedom (virtual nodal quantities), the water pressure and displacement of the virtual nodes are associated with their adjacent nodes through least squares interpolation. Finally, the fully coupled equation in terms of nodal quantities is obtained. The verification cases validate the method. By using this method, the field-scale hydraulic fracturing process is well simulated. The proposed approach is simple and efficient for field-scale hydraulic fracture simulation.  相似文献   

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law          下载免费PDF全文

Ahmad Hajjar  Luc Scholtès  Constantin Oltéan  Michel Antoine Buès 《水文研究》2018,32(16):2510-2525

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.  相似文献   

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials          下载免费PDF全文

Cássio B. Aguiar  Samir Maghous 《国际地质力学数值与分析法杂志》2018,42(16):2018-2046

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.  相似文献   

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock          下载免费PDF全文

N. Zolfaghari  E. Dontsov  A. P. Bunger 《国际地质力学数值与分析法杂志》2018,42(4):587-617

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.  相似文献   

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer     

Evan Earnest  David Boutt  Larry Murdoch  William P. Clement 《水文研究》2019,33(20):2691-2710

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.  相似文献   

不锈钢纤维网对热压铁基预合金胎体抗冲击韧性的影响          下载免费PDF全文

孙吉伟  吴海霞  沈立娜  杨甘生 《探矿工程》2020,47(2):73-77

针对孕镶金刚石钻头在孔底钻进时受到各种冲击力而导致胎体出现裂纹甚至断裂的问题，提出了应用不锈钢纤维网提高胎体抗冲击韧性的方法。通过控制纤维网的摆放位置、目数和层数3个变量，测试了不同组合变量下的胎体抗冲击韧性。通过试验得到：（1）不锈钢纤维网能有效提高胎体的抗冲击韧性，而且X轴位置效果更好；（2）当纤维网的目数相同时，胎体抗冲击韧性随纤维网层数的增加而减小；（3）在纤维网层数相同的情况下，随着纤维网目数的增加，胎体抗冲击韧性呈现先增大后减小的趋势。当添加30目单层纤维网时，胎体的抗冲击韧性达到最大值。通过观察冲击断口、分析胎体断裂机理，解释了出现这种规律的原因，并提出了进一步优化胎体抗冲击韧性的方法。  相似文献   

库车坳陷克深气田致密砂岩储层构造裂缝形成序列与分布规律     

王珂  杨海军  李勇  张荣虎  杨学君  王俊鹏 《大地构造与成矿学》2020,(1):30-46

综合采用岩心、薄片和成像测井等资料,对库车坳陷克深气田白垩系巴什基奇克组构造裂缝的形成序列、分布规律和影响因素进行了研究。结果表明:剪切裂缝和张性裂缝在克深气田均有发育,以直立缝和高角度缝为主,主要形成于近南北向的挤压作用、背斜弯曲拱张作用和异常流体高压作用,微观裂缝切穿胶结物和部分颗粒,早期充填构造裂缝可在后期构造应力和异常流体高压作用下重新裂开成为有效裂缝。克深气田发育3期构造裂缝,其中第3期构造裂缝的形成时间与天然气大量充注期吻合,是工业规模性气藏形成的关键因素。单个断背斜高部位的构造应力低于背斜翼部,因此背斜高部位的构造裂缝线密度相对较低,但背斜弯曲变形使裂缝开度较大,有效性好,单井的无阻流量较高;翼部和断层附近构造裂缝线密度较大,但开度较小,有效性差;构造应力、岩石强度和变形时间的不同造成了构造裂缝特征在各气藏之间具有差异性。储层中部第3砂层组的构造裂缝发育程度中等,充填程度相对较低,并且平面上分布连续,可形成连片分布的储层“甜点”区,应成为克深气田开发中的重点层位。对克深气田构造裂缝起主要贡献作用的为水下分流河道和河口坝微相的粉‒细砂岩。  相似文献   

基于Ruger近似式预测裂缝型储层的可行性研究     

吴姗姗  高刚  桂志先  王鹏  李成毅 《煤田地质与勘探》2020,48(1):189-196

叠前纵波方位各向异性检测裂缝是目前应用最为广泛的方法，该方法主要是利用AVAZ地震资料椭圆拟合的长轴与短轴信息来评价裂缝，该方法在具体应用时，认为裂缝型储层AVAZ响应特征的变化仅由各向异性参数决定，而没有考虑非各向异性参数的影响，从而引起裂缝评价的多解性。为提高该技术的裂缝评价精度，从模型分析入手，将模型中各向异性参数及岩性参数设置为概率密度分布函数，然后采用Monte Carlo随机方法进行叠前AVAZ正演模拟。首先，对各向异性参数对椭圆拟合的影响进行简单分析，得出各向异性参数γ对椭圆扁率B/A与各向异性因子B影响最大、δ次之、ε最小。然后，进一步重点分析不同标准差的速度与密度模型的B/A与B响应特征，得出地层的速度与密度的变化在较小范围内时，该技术评价裂缝才具有一定的可行性，其中，地层纵波和横波速度的变化（尤其是纵波速度）对叠前AVAZ响应影响较大，地层密度几乎没有影响，故当横向岩性变化较大时，该方法在评价裂缝型储层时的结果具有一定的不确定性。   相似文献   

海水海砂混凝土双K断裂参数的确定   总被引：1，自引：1，他引：0  

臧朝会  杨树桐  金亮亮 《海洋工程》2019,37(4):142-150

采用尺寸为515 mm×100 mm×100 mm初始缝高比从0.1变化至0.8的三点弯曲梁试件,利用试验测得起裂荷载Pini、最大荷载Pmax及对应的裂缝口张开位移CMODc等参数,研究了缝高比a0/h对海水海砂混凝土双K断裂参数的影响。并以相同配合比的淡水河砂梁作为对照组,进一步研究了海砂海水混凝土双K断裂参数与普通混凝土之间的关系。结果发现海水海砂混凝土断裂参数在a0/h=0.25～0.7范围内时可以认为是一常数;海水海砂混凝土梁断裂参数比相同配合比情况下的淡水河砂梁大且对边界影响更加敏感。  相似文献